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» ' ABSTRACT 

\o : 

! | The X-ray Pulsar PSR J1811-1925 in the historic supernova remnant Gil. 2-0. 3 has a characteristic age over 
;> 10 times the age of the remnant. This likely implies that its current spin period, 65 ms, is close to its birth 
. spin period. Alternatively, the pulsar may have an unusually high braking index. We report here on regular 
Rossi X-ray Timing Explorer/ Proportional Counting Array {RXTE /PCA) timing observations of the pulsar 
that were designed to measure its braking index. We provide a preliminary phase-coherent timing solution 
which includes a significant v. The braking index we measure is S> 3, likely a result of conventional timing 
\ noise. We also report on a preliminary analysis of the pulsar's unusually hard spectrum: we determine a 
photon index of T = 1.78 ±0.74, for the pulsed component of the spectra, consistent within the uncertainties 
O .■ with previous ASCA and Chandra observations. The pulsed emission of PSR J1811-1925 is seen beyond 
30 keV, and the pulsations remain sinusoidal up to and beyond this energy. 
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INTRODUCTION 

The determination of the ages of pulsars is of crucial importance to neutron star astrophysics. Neutron 
star population synthesis, birthrate, supernova remnant associations, and even cooling timescales (which 
& ', serve as probes of the neutron star equation of state), all depend on accurate determinations of pulsar ages. 
In general, a pulsar's age is estimated by assuming a spin-down evolution of the form 

v oc u n , (1) 

where n is known as the braking index. It can be determined observationally by measuring is, since, as 
is easily shown, n = vvjv 2 . Integrating Equation ^ over the lifetime of the pulsar, assuming a constant 
braking index, results in the following estimate of the pulsar's age: 
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where Pq is the spin period of the pulsar at birth. It is conventional to assume n = 3, appropriate for 
spin-down due to magnetic dipole radiation, and Po << -P to derive the characteristic age r c = Pj (2P). 
Timing noise, stochastic variations in the pulse spin-down (see e.g. Lvne et al.. 19 96). make accurate 



determinations of the braking index impossible for the large majority of the over 1300 radio pulsars 
currently known. The five reliably measured radio pulsar braking indices are in the range 1.4 — 2.9 
| Ka mi ^ al.. 1994bvne et al.. 199ftLvne et al.. IflsdCamilo et al.. 200(1zhang et al.. 200lh . This clearly sig- 



nals a significant difference between pulsar electrodynamics and a si mple dipole spin ning down in vaccuo 



This is of major interest to modelers of pulsar magnetospheres (e.g. iMelatos. 1 997). But it is also very 
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important for neutron star astrophysics in general: if the braking indices of all pulsars are n < 3, then the 
characteristic age systematically underestimates the true age. 

On the other hand, initial spin periods Po that are a non-negligible fraction of the current spin period 
P lead to a systematic over estimation of the true age. The distribution of initial spin periods is also of 
intrinsic scientific interest, since it gives clues to the physics of the formation of neutron stars as well as 
such possibilities as early epochs of rapid spin-down due to gravitational radiation or interactions with fossil 
disks. However, to infer the initial spin period from Equation [21 requires an independent determination of 
the true age, such as an associated historical event or a well determined age from an associated supernova 
remnant. This situation is very rare. The best, and until recently only, example is the Crab pulsar, whose 
measured braking index and known age results in an inferred Po ~ 19 ms. There are very few other historical 
supernovae that might be used. Of course an isolated pulsar's current spin period is an upper limit on the 
initial spin period, so , for example, Po < 16 ms can be trivially inferred for PSR J0537— 6910 in the LMC 
( Marshall et al.. 1998h . 



PSR J1811 1925 The Pulsar at the Heart of Gil. 2 0.3 

Gil. 2-0. 3 is a young supernova remnant plausibly associated with the historical "guest star" witnessed by 
Chi nese astronomers in the year 386 A.D. (|Clark and Stephenson. 1977^ . Radio and X-ray estimates of its 
age ( Creen et al.. 198^Vasisht et al.. 1996h . assuming a type II supernova, support the historical association. 



A young 65-ms pu lsar with spin-dow n energy E = 6.4 x 10 36 erg s 1 was discovered in X-rays by ASCA 
within C1 1.2-0 3 dTorii et al.. 1997h . Recent Chandra observations place the pulsar at the precise center 



of the remnant ([Kaspi et al.. 20011 ). making the association of the two unambiguous. Surprisingly, the P = 
4.4 x 10 -14 , determined from one observation in 199 4 and three observations in 1998, implies a characteristic 
spin-down age of r c = 24,000 yr ( Torii et al.. 1999b . 15 times that inferred from the historical association, 



the remnant radio surface brightness, and the remnant X-ray spectrum. For this large discrepancy to be 
resolved, either the braking index n > 30, or the initial spin period Po ~ 62 ms is very near the current 
spin period. Very r ecently, a similar Po was inferred for the newly discovered pulsar in 3C58/SN1181 
( Murray et al.. 2002h . suggesting that such long initial spin periods may be common. Unfortunately, the 



number of systems where the initial spin period can be estimated is likely to remain small. Therefore, it is 
crucial that all alternative explanations for the age discrepancy observed in Gil. 2-0. 3 be fully explored. 

OBSERVATIONS 

The results presented here were obtained using the Proportional Counter Array (PCA) on board the 
Rossi X-ray Timing Explorer (RXTE). Our observations consist primarily of ~ 20 ks observations taken on 
a monthly basis. 

Phase Coherent Timing 

Deep radio searches have been unable to detect the pulsar ( Crawford et al.. 1998^ , and so X-ray monitoring 



of the source is t he only option for timing. The spin-down rate P of PSR J1811-1925 was first made by 
Torii et al. (1999) using ASCA, essentially on the basis of periods at just two widely spaced epochs. Here 



we present a phase-coherent timing solution for the the pulsar using regular RXTE observations. 

In the timing analysis, we included only those events having energies in the 2-30 keV energy range so as 
to maximize the signal-to-noise ratio of the pulse. Each binned time series was epoch-folded using the best 
estimate frequency determined initially from either a periodogram or Fourier transform (though later folding 
was done using the timing ephemeris determined by maintaining phase coherence; see below). Resulting 
pulse profiles were cross-correlated in the Fourier domain with a high signal-to-noise template created by 
adding phase-aligned profiles from previous observations. The 2-30 keV template is shown in Figure ^ The 
cross-correlation returns an average pulse time-of-arrival (TOA) for each observation corresponding to a 
fixed pulse phase. The pulse phase (ft at any time t can be expressed as a Taylor expansion, 



<f>(t) = <t>(t ) + u Q (t - t ) + \u {t - t ) 2 + -i>o{t - t f + ■■■, 

2 o 



(3) 



where v = 1/P is the pulse frequency, v = dv/dt, etc. and subscript '0' denotes a parameter evaluated at the 
reference epoch t = to. The TOAs were fit to the above polynomial using the pulsar timing software package 
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Parameter 



Value 



v (Hz) 

v (xlO^ 11 Hz s" 1 ) 
v (xHT 20 Hz s" 2 ) 
P (ms) 

P (xl0~ 15 s s" 1 ) 

Epoch (MJD) 

RMS Residual (ms) 

Surface Magnetic Field, B s (xlO 12 G) 

Spin-down Luminosity, E (xlO 36 erg s _1 ) 

Characteristic Age, r c (kyr) 



R.A. (J2000) 

Decl. (J2000) 

Range of MJDs 

No. Pulse Arrival Times 



18 h ll m 29 s .22 
-19° 25' 27."6 
52341-52541 
18 

15.461208393(11) 
-1.00975(28) 
1.40(3) 
64.67799741(44) 
42.240(12) 
52341.8046 



2.7 
1.7 
6.1 

24 



Table 1. Spin parameters for PSR J1811— 1925 derived from phase coherent timing of this source for ~ 9 months. 

TEMPO 1 . Unambiguous pulse numbering is made possible by obtaining monitoring observations spaced so that 
the best-fit model parameters have a small enough uncertainty to allow prediction of the phase of the next 
observation to within ~ 0.2, this was accomplished by two closely spaced observations (within a few of hours 
of each other) followed by one spaced a few days later. We are now capable of maintaining phase coherence 
with regular monitoring observations. 

Arrival time residuals for PSR J1811-1925 are shown in FigureEJ Spin parameters determined via phase- 
coherent timing using nine months of RXTE data are presented in Tabled Phase residuals from the simple 
model given in Table Q are random, though not completely consistent with zero given the uncertainties, 
which are typically ~1 ms. This suggests there might be some pulse "jitter," as is sometimes seen in radio 
pulsar timing. The tabulated spin parameters include a significant frequency second derivative. These 
parameters imply a braking index n = 1980 ± 93. This number is much larger than the canonical value of 3. 
This could reflect highly unusual spin-down behavior, unprecedented in rotation-powered pulsars. If so, the 
initial spin period would not be determined for this pulsar. A more likely scenario is that the large value 
measure for n reflects conventional timing noise, common to rotation-powered pulsars. In this case, a future 
measurement of the true, deterministic value of n is not precluded. 

Spectroscopy 

This source is of particular interest given its relatively hard spectrum, as, unlike most X-ray detected 
pulsars, pulsed emission well above 10 keV is definitely detected. We studied the pulsed spectrum by 
summing all our folded, phase aligned time series into two energy bands (2-10 keV and 10-30 keV). The 
folded profiles in these energy bands were identical. The 2-30 keV profile is displayed in Figure E 

In order to obtain a high signal-to-noise pulsed spectrum for PSR J1811-1925, data from each observing 
epoch were folded at the expected pulse period as was done for the timing analysis. However, here, 8 phase 
bins were used across the pulse. For each phase bin, we maintained a spectral resolution of 128 bins over 
the PCA range. One off-pulse phase bin was used as a background estimator. The pulse profiles were then 
phase aligned, so that the same off-pulse bin was used for background in every case. The remaining phase 
bins were summed, and their spectral bins regrouped into groups of 4 using the FT00L grppha. Energies 
below 2 keV and above 70 keV were ignored, leaving 23 spectral channels for fitting. The regrouped, phase- 
summed data set, along with the background measurement, were used as input to the X-ray spectral fitting 

x http: / /pulsar. princeton.edu / tempo 
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Fig. 1. X-ray pulse profile in the 2-30 keV range. 
Two cycles are shown for clarity. Profiles in the 
2-10 keV and 10-30 keV range are identical. 



Fig. 2. Post-fit residuals for PSR J1811-1925 after 
removal of the best-fit spin parameters in Tabled 



Parameter Value 

f 1.78 ±0.74 

n H (xlO 22 cm" 2 ) 2.36 (fixed) 

Flux (xl0~ 12 erg cm- 2 s" 1 ) 1.18 ± 0.36 (2-10 keV) 

L x (xlO 32 erg s^ 1 ) 2.81 ± 0.86 (2-10 keV) 

Xl V) 0.88 (21) 



Table 2. Spectral fit of the pulsed component of PSR J1811— 1925 to a simple photoelectrical^ absorbed power- 
law model of photon index V. The errors quoted indicate 90% confidence intervals. The X-ray luminosity L x , 
in the 2—30 keV range, was calculated assuming a distance of 5 kpc (Becker et al., 1985; Green, 1998) and 
assuming a beaming fraction of 1 sr. 

software package XSPEC 2 . Response matrices were created using the FTOOLs xtef ilt and pcarsp. We fit 
the data using a sim ple photoelectrically absorbed power law, holding only Nh fixed at the value found by 
Roberts et al. (20021 ) (see Table |2J). Uncertainties were measured using the XSPEC command steppar. The 
photon index T we o btained is consistent within the uncertainties with the value (T = I.II^q'^) found by 
iRoberts et al. (20021 ) using Chandra observations. 

DISCUSSION 

The ratio of the pulsed X-ray luminosity in the 2-10 keV band to E gives an efficiency, 

r? x = h = 0.06%, (4) 



assuming isotropic emission. Comparing the efficiencies calculated by iPossenti et al. (2002) for the 41 



known rotation-powe red X-ray pulsars, PSR J1811-1925 is the third most efficient , with PSR J1846-0258 
( Gotthelf et al.. 20 00") being the first, and PSR B 1509-58 dCusumano et al.. 200lh the second. Note that 



the efficiencies calculated hv TlWnti et ah ttnoft often included unpulsed contributions resulting from pul 



sar wind nebular emission, therefore their efficiencies can only be considered upper limits to the pulsed 
efficiencies. It is interesting to note the similarities among the three most X-ray efficient rotation-powered 

2 http://xspec. gsfc.nasa.gov 
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Fig. 3. Top: Phase-averaged spectrum of the pulsed component of PSR J1811— 1925. The dash-line indicates a 
simple photoelectrically absorbed power-law model. Bottom: residuals after subtraction of the model. 



pulsars: (i) they all have long duty-cycles, with profiles that remain broad even to high energies (> 10 keV) 
(ii) they have very hard spectra with photon indices T ~ 1.1 — 1.3 (iii) they are all very young (< 2000 yrs) 
pulsars in the centers of supernova remnants (iv) they all have pulsar wind nebulae. Furthermore, none of 
these pulsars have been detected in hig h energy T-rays (E > 100 MeV). This is in contrast to pulsars such 



as the Crab, Vela, and PSR B1706-44 ( Gotthelf et al.. 20021 ') . which have very narrow X-ray pulse profiles, 



relatively low X-ray efficiencies, but are all detected in high energy 7-rays. 

The spectrum and pulse morphology can give us a clue to the emission mechanism at work in PSR J1811- 
1925. High-energy emission mechanism models generally fall into two classes: polar cap models (e.g. 
Harding et al.. 2 002) and outer gap models (e.g. iRomani. 20021 ). The two harder components of the four 



component X-ray profile (in the RXTE band) of the Vela pulsar, which are as hard or h arder that PSR J1811 - 
1925 are due to particle curvature emission according to the polar cap model (see H arding et al.. 2002h . 



These two components are very narrow and phase aligned with the double peaked 7-ray pulse. For 
PSR J181 1-1925, the broad profile and the absence of 7-ray pulsations make a similar mechanism unlikely. 
In the outer gap model, t he broad, hard p rofile of PSR J1811-1925 would be from the tail of the synchrotron 



Ipi 

component. According to R omani (20021 ). the peak of the synchrotron spectrum for PSR J1811-1925 would 



occur at E ~ 31 MeV, this could be verified with 7-ray observatories such as GLAST and INTEGRAL. 
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